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Abstract Recent experiments 1 on aligned DNA show hexatic order with 
no sign of macroscopic chirality. I make the analogy between smectic 
liquid crystals and chiral hexatics and show how the absence of chirality 
cannot occur in a thermodynamic phase of chiral molecules. In addition, 
I discuss the microscopic origin of chiral mesophases in liquid crystals 
and show that, within the context of central forces between "atoms" on 
"molecules" , chiral interactions can occur only if there are biaxial corre- 
lations between the mesogens. Weak biaxial correlations can therefore 
lead to small cholesteric pitches. 

(5 June 1996) 

INTRODUCTION 

Experimental realizations of the liquid crystalline hexatic, a phase with bro- 
ken orientational but not translational order 2 , are difficult to find 3 . More exotic 
hexatic structures have been proposed 4 in three dimensions in which there is ne- 
matic order and, simultaneously, hexatic order in the plane perpendicular to the 
director n. Until recently, only the liquid state (the regular nematic) and the 
hexagonal state (the hexagonal columnar phase) had been observed, typically in 
long, chiral molecules such as DNA 5 . In fact, the liquid-like phase of chiral poly- 
mers is cholesteric - the nematic phase will always start to twist in the presence 
of chirality. The hexagonal columnar phase, however, need not twist: there is a 
thermodynamic region of elastic constants in which twist is expelled. The twist 
can come in either via a tilt-grain-boundary state in which the nematic order twists 
or via a moire state in which the polymers braid about each other and the local 
crystalline axes rotate along the polymer axis 6 . 

In a new set of experiments 1 a phase with both nematic and hexatic order was 
found in DNA, in addition to a cholesteric phase at lower concentration. In the first 
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section of this talk I will argue that this is rather surprising: in an "N+6" phase 
Landau theory predicts that, unless a microscopic parameter is precisely tuned, 
either the nematic order or the bond-orientational order must twist. If that were 
the case, the X-ray scattering in the plane perpendicular to the nematic director 
would be a powder average over many different, rotated hexatic regions. Thus one 
might expect that there should be a ring in the q ± plane, rather than the observed 
cos66> modulation as shown in Figure 1. 




I 

FIGURE 1. X-ray structure function in the plane per- 
pendicular to the nematic direction 1 . This diffraction 
pattern contains a non-zero cos 69 component and no 
measurable cos6n# for n > 2. The small amount of 
cos 29 can be attributed to the misalignment of the X- 
ray beam. (Figure provided courtesy of R. Podgornik). 

If the chiral interactions are small, however, the hexatic order may not twist. 
Fortunately, there are reasons to believe that this may be the case: in the last sec- 
tion of this talk I will argue that molecular chirality manifests itself as mesoscopic 
chirality (i.e., cholesteric pitch) via a subtle short-range correlation of molecular 
orientation. Naively, one might expect that the cholesteric pitch should be on the 
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order of the molecular scale divided by a fraction of a radian, P 20A/(0.l7r). 
Typically pitches are much longer, on the order of microns. The necessity of in- 
tramolecular correlations may shed some light on the "naturalness" of observed 
cholesteric structures. 

CHIRAL BOND ORDER AND EXPULSION OF CHQLESTERIC TWIST 

Based on symmetry let us construct a Landau free energy to describe the 
phases of chiral molecules with nematic and hexatic phases 6 ' 7 . In the non-chiral 
nematic phase, orientational fluctuations are controlled via the Frank free energy 
density: 

T n = ^ (V-n) 2 + *± (n-V x n) 2 + ^ [n x (V x n)] 2 . (1) 

If there is hexatic order, there will also be a spin stiffness for the hexatic bond- 
angle 6q. When determining the free energy for 6*6, one must take into account 
the transformation properties of 9q under the nematic symmetry n — > — n. Since 
angular changes must be measured with respect to some vector (via a "right-hand" 
rule), under nematic inversion 9q — > — 9q. This does not change the sense of any 
twisting present in 9q since the angular change is measured with respect to n. 
Taking into account the local nematic anisotropy, the elasticity for 9 6 is: 

^ = ^ T ^[n-V^] 2 + ^[V^] 2 (2) 
+ K 6 V x n • V#6 + K'q (n-V x n) (n- V6> 6 ) . 

These terms are invariant under the simultaneous change (n, 9q) — > — (n, 9q). The 
combined free energies have been considered before and constitute the elasticity 
theory for the "N+6" phase 4 ' 8 . 

In a system composed of chiral molecules (such as DNA), one expects addi- 
tional chiral terms invariant under nematic inversion but which change sign under 
spatial inversion. Since nematic symmetry already forces all terms to have even 
powers of n and 9q, a character need not be assigned to these variables (i.e. vector 
or pseudovector, pseudoscalar or scalar, respectively). Thus I include the following 
chiral terms, quadratic in the fields: 

r = K 2 q n-V x n + K%h-V9 6 . (3) 
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The first term is the usual chiral term which favors a cholesteric texture n = 
[cosqoz, singo-2, 0] while the second term favors the twisting of the hexatic order 7 . 
The second term has been considered before in cholesteric melts 9 and chiral polymer 
crystals 6 . Similar theories could be constructed for molecules or molecular packings 
with other sorts of in- plane symmetries, for instance tetrahedral molecules 10 . 



FIGURE 2. A chiral hexatic. In each plane the bond- 
angle order parameter is 6q = 0$ mod 27r/6. The bond- 
angle rotates smoothly between the planes. Each plane 
of constant bond-angle is analogous to a smectic layer. 
This phase has no density modulations 



In the following I simplify to the case = K A = Ka, Kq = K' e = and 
employ the ansatz first considered by Meyer 11 for a nematic in a magnetic field: 
n = [cos <p cos qz, cos <p sin qz, sin (p], with (f> constant. Minimizing the total free 
energy F = J d 3 x \T n + T% + T*] with respect to 9q yields an effective free energy 
for n (note that the mean-field equations only allow constant V#6) : 





(4) 
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This is precisely the free energy considered by Meyer for a nematic with negative 
magnetic susceptibility in a magnetic field. It is straightforward to show that for 
K 2 < K 3 the cholesteric phase (0 = 0) is the lowest energy state (within this 
class of solutions 12 when K 2 q^ > Ka^q and that the nematic phase (0 = ir/2) is 
lowest in energy otherwise. When K 2 > K 3 an intermediate conical phase with 
< sin^ < 1 is allowed 11 ' 7 , though typically this inequality of the Frank constants 
is not satisfied. 

We have just seen that the presence of hexatic order can act as a field. In fact, 
the chiral coupling to the hexatic bond-order field makes this system identical to 
that of a smectic liquid crystal composed of chiral mesogens with a layer spacing 
of 27r/(6(/o)- The hexatic "wave" plays the role of the smectic layers. Each "layer" 
is a place where the bond-angle is 9$ = mod 2tt/6. This is illustrated in Figure 2. 
As in smectics, there will be a myriad of phases, though, in analogy with a type-II 
superconductor 13 , the simplest phases will be the cholesteric phase which has no 
chiral bond order (normal metal), the chiral hexatic phase with no cholesteric twist 
(Meissner phase), and finally the twist-grain-boundary (TGB) phase 14 (Abrikosov 
flux lattice phase). 

This analogy can be made precise 7 by introducing the hexatic order parameter 
ifjQ. The Landau theory which describes the liquid-to-hexatic transition is: 

^6 = |W> 6 | 2 + r#6| 2 + ^6| 4 , (5) 

where r oc (T — T c ). In the ordered phase (T < Tq) ipe = I ^6 1 e" 16 * 6 and \iPg\ 2 = 
Ka/72. Adding the second term of (3) as well as the Frank free energy for a 
nematic with pitch q (^[n; q ]) yields: 

Tl = |W 6 | 2 + r# 6 | 2 + t# 6 | 4 - 6iq h-(ip*ViP - W") + ?o]- (6) 

Completing the square, this leads to: 

Tl = | (V + 6zQ n) VI 2 + (r - 36<? 2 ) |V>| 2 + «|V>| 4 + H™, ?o], ( 7 ) 

which is precisely the free energy of a smectic-A liquid crystal with layer spacing 13 
a = 2n/(6q ). 

Drawing upon the smectic-A analogy, one would expect that cholesteric twist 
will be excluded in favor of chiral hexatic order. In particular, this means that from 
plane to plane the sixfold order should rotate if there is nematic alignment. This 
is not what is seen experimentally - there are six distinct spots in the scattering 
perpendicular to the nematic director 1 . While the molecular symmetry suggests 
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that the chiral terms in (3) must be present, they may nonetheless be small. In the 
next section I will analyze the origin of chiral interactions among molecules and 
argue that they could be anomalously small. It should be pointed out that if all 
the chiral terms were set to zero there is no a priori reason why the phase sequence 
nematic-N+6-hexagonal-columnar-phase could not exist 15 . 

THE MICROSCOPIC BASIS OF CHIRAL INTERACTIONS 

A molecule is chiral if its symmetry group does not contain the element S n : a 
rotation around a C n -axis by 27r/n followed by a mirror through the perpendicular 
plane 16 . Moreover, all chiral molecules must have n-axial order around some axis 
since the groups and are not allowed for any real objects 17 : the only infinite 
subgroups of the rotation group in three dimensions are Coo V and Pqo/i, both of 
which contain Si. As we shall see, the lack of uniaxial symmetry plays a crucial 
role in the chiral interactions. Let us first concentrate on calculating the cholesteric 
pitch qo, or, more accurately, K 2 qo the coefficient of (n • V x n), and return later 
to estimating %. The following will only be a sketch of the main results presented 
elsewhere 18 . To this end, let us consider two planes of molecules separated along 
the x-axis by R, the first aligned along the z axis, the second along (z cos 9+y sin 9) . 
Expanding for 9 <C 1 yields: 

F = j d 3 x K 2 q djn y + 0{9 2 ) « K 2 q L y L z 9, (8) 

where L y and L z are the sample dimensions. There is a famous expression for the 
coefficient of 9 in this expansion-the torque. Thus L y L z K2qo is the torque that one 
plane of molecules exerts on the next plane of molecules when they are aligned along 
a common axis. In the simplest case, let us consider weakly interacting molecules 
and thus build the torque up out of nearest neighbor, pairwise interactions across 
the planes. If there are N molecules per plane then the torque between two of the 
molecules is AK 2 q , where A is the free area of each molecule. 

We see that the problem of calculating K 2 q reduces to the problem of calcu- 
lating the torque between two molecules separated by a distance R when they are 
aligned with each other and perpendicular to their separation. For simplicity, let us 
assume that the atoms on the molecules all interact via the same central potential. 
This analysis necessarily excludes quantum interactions which can be non-local 
and which do produce chiral terms in the Landau theory of liquid crystals 19 , but 
does include, as a limiting case, hard-core steric interactions which can be built up 
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out of short range central forces. The relevant component of the torque is along 
R = R/|R|, where R is the separation between the centers of mass of each molecule 
and R ■ n = 0. The atoms on each molecule are at a position r a relative to the 
center of mass and a runs over the atoms on the molecule. Assuming a central 
potential V(R), the relevant component of the torque is: 

* r = -jfjj E R l ^jkr J p d k V(R + rp- r Q ), (9) 

' ' a/3 

where a runs over the atoms on the first molecule and (3 runs over the atoms on 
the second molecule. An expansion in powers of (|r|/|R|) is certainly reasonable in 




FIGURE 3. Prototypical chiral molecule. The spheres 
indicate the positions of the atoms. When the angle 
7 is mr/2 this molecule is not chiral. The chiral order 
parameter for this molecule is ip = — (l/8)w 4 L sin(27), 
where L and w are the dimensions shown and L > w 
(not drawn to scale). 

the dilute limit and the general symmetry of the interaction should not be signifi- 
cantly altered at higher densities. Note that terms in the expansion which are even 
in r must vanish: under r — > — r both molecules will change their handedness and 



7 



thus the chiral pitch will change sign. Note that d k V(R) oc R k and thus from the 
antisymmetry of the first order term will vanish. The third order term is: 

(R • r) = -~ (rj? - C) (rj, - r< ) d fc d^(R). (10) 

' ' 'a/3 

The potential term d k d m d l V(R) will, by symmetry, have only two types of terms: 
those proportional to 1) R k R m R l and 2) (R k S ml + R m S kl + R l S km ). Again, due 
to the antisymmetry of e^fc, terms of the first type will vanish and the only part 
of (10) which survives is: 

(* ' r ) 3 « E ■ (t p - r«) . (11) 

a/3 

Note that every term in (11) will contain the single sum which is identically 

since the r M are measured from the center of mass. Thus, the third order term 
will vanish. 

Fortunately, perseverance will pay off and the fifth order term will not identi- 
cally vanish. At fifth order there are a variety of terms. Again, any term which is 
a sum over a single factor of will vanish. This leaves only terms with three r % 
on one molecule and two r 7 on the other. Taking the average yields 18 : 

(R • t) = U(R)Q u e ljk {s 3 2 lm B km + S{ lm B km } , (12) 

where U(R) = {\R?f^{R 2 /2) - \Rf^\R 2 /2)} is a function of R = |R|, V(R) = 
/(.R 2 /2), Q = \n l ri — (l/3)S l ] is the nematic alignment tensor (equal for both 
molecules), the two tensors S^ lm and B km are defined on molecule p = 1, 2 via 

^ m = e f^i - ^m 2 i^j + + } 

A*6p ^ - J 

and 5 lJ = [5 % i — n*n J ] is the transverse projection operator. The tensor B km 
measures the biaxial orientation of the molecules in the plane perpendicular to the 
nematic direction, which is assumed to be the principal axis of the molecule with 
the largest moment. The other tensor 5 ,j7m is some measure of the chirality of 
the molecule, but is not a measure by itself: indeed, some components of 5 ,j7m 
will be non-zero for achiral molecules and zero for chiral molecules. Consider, 
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however, the order parameter ip p = S J p lrn eijkQ %l B km on a single molecule p: in 
the coordinate system of the principal axes, ifj = (X y — X x ) ^^r^r^r^ where A x 
and A y are the eigenvalues along the x and y directions of the single-molecule 
moment of inertia tensor and z is the nematic axis. This order parameter vanishes 
identically for achiral molecules, since under inversion the principal axes do not 
change but the sum will change sign. Solving for S p lm = ifj p €ijkQ tl B km / (2B 2 ) + 
(symmetric permutations), (12) reduces to: 

(K.r) 5 = ^U(R)^^. (14) 

where B 2 = B\ = B\ is the same on the two identical molecules. Up to this point, 
the average over the n-axial direction has not been performed. In any mean field 
theory ( B km ) will either be zero (uniaxial) or non-zero. In the latter case the 
entire phase will be biaxial, with long-range biaxial order. If the same mesogens 
in racemic mixtures have only uniaxial nematic phases then (B km ) = (i.e., 
no long-range biaxial order). Thus if the molecules rotate independently around 
their nematic axes then the average ( Sf m i?| m ) = ( B km ) ( B krn ) = 0. Thus in a 
uniaxial phase, chiral interactions will come about only via short-range correlations 
of biaxial order. In other words, upon averaging, (14) becomes 

/ \ / flkm ~n>km \ 

( R • r ) 5 = V U(R) K - 1 B2 2 1 = V U(R)g B (R), (15) 

where Qb{R) ~ e~ R ^ B is the normalized biaxial correlation function, and £_g is 
the biaxial correlation length, which one might naively expect to be on the scale 
of the intramolecular spacing. 

The previous discussion shows that there are two factors which lead to a large 
chiral parameter qo: the strength of the chiral order parameter ip on the molecule 
and the range of the biaxial correlations. If either is small, it is reasonable to 
expect a very long pitch. Note that almost cylindrical, corkscrew-like molecules 
are not very biaxial and thus ip would be small. In addition, in the DNA system 
studied at NIH, the molecules are on the order of 40A apart (center-to-center) 
and are each roughly 20A wide. Moreover, the Debye screening length at the 
typical salt concentrations investigated is on the order of 5A and thus one might 
expect that the biaxial interactions would be weak and hence £s to be small. 
Together, the smallness of both factors suggests that the chirality might be very 
small. In fact, it is even possible that in the hexatic phase biaxial correlations 
are smaller than in the cholesteric phase, despite the higher molecular density. 
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Since the cholesteric phase picks out a preferred, biaxial direction (the pitch axis) 
the biaxial correlations could be longer-ranged (even possibly long-range). In a 
hexatic phase, however, the biaxial order parameter might be suppressed due to 
the 6-fold "crystal" field of the hexatic. Note that if the biaxial correlations are 
removed there could still be "6-axial" correlations. In that case, however, the 
analogous order parameter to ip would be more complicated and would include 
higher geometric moments. Moreover, the expansion of R-r in powers of (|r|/|R|) 
would vanish, upon averaging, even at fifth order, and one would have to go to, at 
least, ninth order which would lead to a very small interaction at dilute separations. 

Finally, let us estimate the amount of twisting of the hexatic order by consid- 
ering a simplified two-polymer model. Consider one polymer helically winding its 
way around another straight polymer at a mean separation a: the rate at which the 
separation vector rotates along the straight polymer axis is a simple estimate of % 
and will be precisely equal to the pitch of the helical polymer trajectory. There are 
two major contributions to the energy of the helical polymer: the first is the usual 
bending energy of a stiff polymer, while the second is the "nematic" free energy 
for the two molecules which includes the cholesteric term proportional to go- Take 
the first polymer to follow the curve r(z) = [acosqoz, asinq^z, z] and the second 
to lie along the 5-axis at x = y = 0. In this case n • V x n « (aq )/a = and so 
the free energy of the second polymer is: 

F = J dz j| {1 + (ago) 2 }" 372 a 2 (g~ ) 4 + ^ [q + qo?] • (16) 

Minimizing for small ago and taking the bending stiffness k = k B TLp, where Lp 
is the persistence length, yields: 

3 K 2 K 2 
2L P (g ) + — g = —j^Qo- ( 17 ) 

Taking the typical DNA values for K 2 = 10~ 6 dyne, g = 2n / 1 jjam and L P = 600A, 
I expect that at room temperature go ~ —go- It should be noted, however, that as 
Lp — > oo, go — > 0, in accord with one's intuition - as the polymers become infinitely 
stiff there can be no hexatic twisting. More generally, the elastic constant K 2 has 
been estimated for semi-flexible polymers with mean spacing a (roughly 40A in the 
DNA samples) to be 20 K 2 w k B TLp /S /a 4 / 3 . Substituting into (17) I find g g 

1 /3 

for pitches P = 2n/q > 2n (Lpa 2 ) w 600A, which is almost always the case. 
Thus, because the typical pitches are much longer than any intrinsic molecular 
length, the polymer stiffness plays no significant role in reducing g . 
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